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Classical and quantum conduction on a bond-diluted Bethe lattice is considered. The bond 
dilution is subject to the constraint that every occupied bond must have at least fc — 1 neighboring 
occupied bonds, i.e. fc-core diluted. In the classical case, we find the onset of conduction for fc = 2 
is continuous, while for fc = 3, the onset of conduction is discontinuous with the geometric random 
first-order phase transition driving the conduction transition. In the quantum case, treating each 
occupied bond as a random scatterer, we find for fc = 3 that the random first-order phase transition 
in the geometry also drives the onset of quantum conduction giving rise to a new universahty class 
of Anderson localization transitions. 



I. INTRODUCTION 



Theoretical study of the interplay between uncorrelated disorder and quantum mechanics began with the seminal 
work of P. W. Anderson more than fifty years ago [l| . Anderson predicted the existence of spatially localized single 
particle states provided there is sufficient disorder in the potential. This finding identified the possibility of a phase 
C . transition from conducting (extended states) to insulating (localized states) with increasing variance in the distribution 
'T^ ' of disorder. The nature of the metal-insulator transition eventually became more transparent with the introduction 
C^ , of the phenomenological one-parameter scaling theory of localization [2] based on ideas developed by Wegner Q and 
^ Thouless [4|. The scaling theory was subsequently fortified by a nonlinear sigma field theory and a perturbation 
theory near two-dimensions allowing for calculation of the set of exponents characterizing the continuous Anderson 
Cd ' localization transition 0-0] • 

There exists another approach to studying the interplay between quantum mechanics and disorder via the model 

T 1 of quantum site percolation 8] . Consider a binary alloy where the potential landscape is represented by two different 

' ^ ] energies randomly distributed throughout the system with probabilities p and 1 — p respectively. In the limit where 

the energy difference approaches infinity, the quantum particle can only access (hop to) one of the two types of atoms. 

These randomly distributed accessible sites may or not span the system and, hence, affect the quantum conduction. 

If the accessible sites do not span the system then surely the system is insulating, for example. 

Some controversy surrounds the study of quantum percolation. For instance, some argue that quantum percolation 
transition in the same universality class as the Anderson transition [9|, [lOJ. However, others find evidence for a 
transition in two dimensions [llj^lSi], contrary to the scaling theory. One issue that appears to be resolved to some 
h — [ ' extent in three dimensions is the fact that Pq > Pc, where Pc signals the onset of the geometric percolation transition 
\^ . and Pq signals the onset of extended single particle wavefunctions '14] . 

t:]" ' The Anderson model and quantum percolation have been analyzed on the Bethe lattice |15l - [l8| . The loopless 

ly-^ structure of the Bethe lattice makes it amenable for analytic study, hence, it will be implemented here. For instance, 
f~^ ' Abou-Chacra and collaborators were able to obtain closed form expression for the breakdown of localized states 
f^ I in terms of the value of the potential, the width of the disorder and the coordination number of the lattice [15| . 
T-H . Harris |17l [l8| analyzed a bond version of quantum percolation and found that the exponent associated with the 
divergence of the average finite cluster size as the transition is approached from below is the same as classical 
percolation but with pq > Pc- 

^\ ' Traditionally, the Anderson model and quantum percolation are models of quantum transport with short-ranged, 

'^ , uncorrelated disorder. More recently, eigenfunction studies of power-law diluted chains have been conducted 19]. 
More specifically, the hopping probability of an occupied bond scales with, r, the distance along the chain, as l/r^+'". 
As V is decreased below 0.68, extended states emerge presumably due to the fact that the system is becoming mean- 
field-like such that this result is not contrary to the scaling theory of localization. 

How is the Anderson transition or the quantum percolation transition affected by other types of disorder — for 
example, correlated disorder where the correlations arise via local constraints on the occupation of bonds? The 
simplest model of correlated percolation is fc-core bond percolation where every occupied bond must have at least 
fc — 1 occupied neighboring bonds |20l - |24| . To enforce this constraint, bonds are initially occupied independently and 
at random with probability p. Then, those occupied bonds with less than fc — 1 occupied neighboring bonds are 
rendered unoccupied. This removal procedure proceeds recursively throughout the lattice until all occupied bonds 
satisfy the fc-core constraint. Please see Fig. [T] for an example on the Bethe lattice. 

As for the fc-core percolation transition, in mean field, for fc > 3, the fraction of occupied bonds in the spanning 
cluster. Poo, is finite at the transition. This result is to be contrasted with fc < 3, where the fraction of occupied 



bonds is zero at the transition. While the fc > 3 transition is discontinuous in terms of the order parameter, Poo, there 
exists at least two diverging lengthscales exhibiting evidence of a random first-order phase transition Q^ . Therefore, 
the A; > 3 represents a new universality class differing from the ordinary, uncorrelated percolation model. 

Given this new universality class in the geometric percolation transition from disconnected to connected due to 
constraints on the disorder, let us return to the theme of the interplay between disorder and quantum mechanics. 
How does the random first-order phase transition in the geometry affect the onset of quantum conduction? Could 
the discontinuous nature of the nature allow for pc = Pq as well as allow for the discontinuous onset of conduction, 
providing evidence for a new universality class in quantum localization transitions? 

We will provide an answer to this question via analysis of quantum conduction on a fc-core diluted Bethe lattice. 
Before doing so, we present results of classical conduction on a k-core diluted Bethe lattice since we will implement some 
the machinery in the quantum limit as well. Note that the k — 1 case has been analyzed previously by Stinchcombe [25| 
and Kogut [26[ has analyzed a site version of the k = 3 case. Therefore, the paper is organized as follows: Section II 
provides the classical analysis. Section III provides the quantum analysis using Landauer conduction, and Section IV 
discusses the implications of our results. 




FIG. 1: Here fc = 3 and z — 4. The shaded circles denote branches that are fc — 1 connected to infinity. The removal of bonds 
1 and 2 eventually triggers removal of bond 3 and bonds emanating from vertex a, including the shaded circles. The remaining 
three branches emanating from the center site survive the removal process. 



II. CLASSICAL fc-CORE CONDUCTION 



A. Geometry of fc-core clusters 



Consider a seed vertex from which a lattice with coordination number z and N generations emerges. Each generation 
is constructed by recursively adding z — 1 bonds to an mth generation site forbidding the formation of loops to produce 
a rooted Bethe lattice. Consider the missing z bond of the seed site to survive the fc-core removal process as specified 
in the introduction. Then the entire connected cluster of occupied bonds, each occupied with probability p, will 
survive the fc-core removal procedure if each occupied bond in the mth generation has fc — 1 occupied neighbors bonds 
in the {m + l)th generation. To determine the percolation properties of such geometrical structures, we define R, the 
probability that an arbitrarily chosen branch leaving a given site is not in an infinite fc-core cluster with 



R=l -p+p 



fe-2 

E 

n=0 



1 



i?^ 



-"(1-i?)". 



(1) 



The arbitrarily chosen branch is not part of an infinite fc-core cluster if (a) the bond is not occupied or (b) the bond 
is occupied, but less than fc — 1 of its neighboring bonds are occupied. For k = 2, the equation for R reduces to the 
ordinary bond percolation problem with R decreasing continuously from unity just above the transition. In particular, 
for z = 3 with p — pc = e << 1, 

( 1 if e < 

^ - i ^ - 1 = t4- - 1 ~ 1 - 4e if e > 0. ^2j 

For fc > 3, however, R jumps discontinuously from unity at the transition such that fc > 3 represents a different 
universality class from ordinary percolation. For example, for fc = 3, z = 4, 

1 if e < 

lot 

where Pc = |, i?c == j, and Rq = ^. 



^ " ^ R,- Roe^/^ if e > ^^^ 



B. Classical conduction formulas 

Each occupied bond denotes a conductor with conductivity ctq > 0, while each unoccupied bond denotes a conductor 
with zero conductivity. The probability of an arbitrarily chosen branch leaving a given site having a conductivity a 
is given by 

P{a)^R6{a) + {l-R)H{<j), (4) 

where 

H{a)^pJ2 (^^^){l-RrR''''" fdai...danH{a,)...H{an)6{a-Sn), (5) 



with 



Co + Tn 

and 



(6) 



r„ = ^a.. (7) 



i=l 

Note that Sn assumes that the bond starting off the arbitrarily chosen branch is occupied. We have also invoked the 
self-similarity of the Bethc lattice such that H^a) is equivalent from one generation to the next. 

Once H(a) is determined for an arbitrarily chosen branch, the microscopic conductivity distribution for the system 
can be computed via 

p('t) = ^ Q (1 - RrR'-"5{a) + J2 P) (1 - RTR'-" f d^i- da„i/(ai)... if (a„)(5(a - r„). (8) 

n=0 ^ ^ n=k ^ ^ "^ 

The first term represents those realizations where the fc-core criterion is not met for the z branches and the second for 
those realizations otherwise. We note that P{a) is normalized as is H{a) and p{a). Finally, the average microscopic 
conductivity of system, < a >, is given by 

< (T >^ daap{<T). (9) 



C. k = 2,z — 3 case 

We first compute H{a) determined by 

H{cr) =2pR f daiH{ai)S{(T - Si) + p{l - R) ( daid<72H{ai)H{a2)5{<7 - S2) (10) 

with R given by Eq.[21 One expects the onset of nonzero average conduction to occur at p^i the occupation probability 
above which there exists a spanning cluster with probability unity. Since the geometric transition is continuous, we 
expect that the conduction transition is also continuous. Therefore, we propose the ansatz, 

H{a) ^ ^H{^). (11) 

We set 9 — 1. Inserting this ansatz into Eg. [TOl yields 

H{a) = 2(. + i)(l - 4.) / daiHC-l)Sia - ^^) + (e + \)Ae f daida2Hi^)H{^)S{'r i^l±^). (12) 
2 J e CTi + ctq 2 J e e cri + 0-2 + ctq 

We let - = a; and — — Xi for i = 1, 2, multiply both sides by exp(— sx), and integrate over x from zero to infinity to 
arrive at 

H(s) = 2(e+i)(l-46) fdxiH{xi)eM-T^) + i^ + h^e j dxidx2H{xi)H{x2)eM- , "^f +";^^ ), (13) 

^ J "•" ~ ^ J ^+ -^\Xl+X2) 

where 'H(s) is the Laplace transform of H{x). 

Expanding in powers of e and collecting terms, the zeroth term in e yields the identity 



'H(s) = / dxiiJ(a;i)exp(-sa;i). (14) 

The first order term yields 

= — 2 / (ixii/(xi) exp(— sxi) H / (ia;ia;i_ff(xi) exp(— sxi) + 2 / dxidx2H{xi)H{x2)cy~'p{—s{xi+X2)) (15) 

J CTQ J J 



leading to the differential equation. 



H"{s) = 2^n(s){l-n{s)). (16) 

s 



As for the two boundary conditions, the normalization on H translates to 'H(O) = 1. We also choose H{oo) — 0. 
We solve this system of equations numerically using Mathematica [27:]. A result is shown in Figure [5] The large s 
behavior is determined by 'W'is) — 2^'H{s) with 'H(s) ^ exp(— 2-v/2/(7n5^^^)- We note that the same form of the 
differential equation emerged from the analysis by Stinchcombe [25| to second order in e. (See also similar analysis 



by Heinrichs and Kumar [28 1. 



The above analysis justifies the ansatz for H{a) for e << 1. Inserting the ansatz into the expression for < ct > to 
determine the scaling of < u > with e yields 

<a>=96 f dxxH{x)e^ +0{e^). (17) 

This result is to be contrasted with the result for fc = 1 analyzed by Stinchcombe (25| where quadratic scaling with 
e was computed. This difference is due to the fact that while the expression for R is equivalent for k — 1 and k — 2, 
the fraction of bonds participating in the infinite cluster. Poo, scales differently with e for A; = 1 and k = 2. For fe = 1, 
Poo ^ e, while for k — 2, Poo ^ e^ such that the conductivity should scale with e as above. This scaling holds for larger 
z as well. We note that the fc = 2 case is equivalent to biconnected percolation studied by Harris ^9]. Harris found 
that the exponent for ji/i-connectedness is given by jf3, where /3 = 1 is the order parameter exponent for ordinary 
percolation. Presumably, the conductivity exponent generalizes to jf3 + 1 in the ji/i-connected case, which is different 
from fc-core in that there is no culling process. 




FIG. 2: -H(s) for ctq = 1/2. 



D. k = 3, z = 4: case 



The equation for H(a) now reads 

H{a) = 3pRil - R) f dc7ida2H{ai)H{a2)Sia - S2) + ?{! - Rf j daida2d(j^H{ai)H{a2)H{c7:i)5{(j - S^a), (18) 

where R is now given by Eq. [3l Since there is a jump in R at the transition for this value of fc, one expects the 
transition in the conductivity to be discontinuous as welL We propose the foUowing scahng form for H{a): 



H{a) 







if e<0 



Hc{cr) + e^K{a) if e > 0. 



(19) 



The normalization on H{a) implies / daHdcr) ~ 1 and J daK(a) = 0. 
Using Eq. \W\ and expanding in e, we arrive at 

Hci<y) = ^{f daidc72H,iai)H,{a2)Sia - S2) + J daida2da3H,{ai)H,ia2)H,{a3)S{a - S3)) (20) 

for the e independent terms. If A > 1/2, then the terms of order e^/^ yield 

H,{a) = j daida2H,{ai)H,{a2)5{a - S2) (21) 



Hc{(t) = / d(Tld(T2d(T3Hc{(Tl)Hc{(72)Hc{(T3)S{a ~ S3) 



(22) 



Equation [5T] yields solutions, Hc{ct) — 0,S{a), or S{a — ^). The last solution is the only one allowed given the 
construction of Hc{a), however, this solution conflicts with the solution Hc{cr) ~ 6{(t — Ictq) from Eq. [221 So, A < i. 
If A < ^ , the terms of order e'*' imply 



K{a) - f daida2H,{a,)K{a2)S{a -82) + ^ j da^da2da3H,{a^)H,{a2)K{a3)5{a - ^3). 



(23) 



However, the linear integral equation for K{a) dictates that K{a) can be arbitarily rescaled by a factor. Since any 
physical solution should be unique, K{<t) == is the only solution. Therefore, we rule out A < i. This leaves A = i. 
Then the terms of order e^/^ are 



Kia) - I da^d<j2H,{ai)[K{a2) - {'iV2/^)H,{(j2)]5{(j - S2) 
[daida2da3H,{ai)H,ia2)[isV2/4:)H,{a3) + (3/2)if (a3)](5(a - ^3). 



(24) 
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FIG. 3: Left: Plot of Hdp-) with binsize 0.0025 for k = 3, z = 4 classical bond percolation. Right: Plot of K{a) for the same 
conditions. Hc{(t) converges after 14 iterations and K(o) after 6 iterations. 



To justify the above ansatz, we must find nontrivial solutions for Eqs.[20land [23] with A = i. We do so numerically 
by making an initial guess for Hc{(y) and K{a) and solving the integral equations iteratively until both equations are 
satisfied within some tolerance. This algorithm is first implemented for Hc{<j) and then that numerical solution is 
used to solve for K{a). For the rest of the classical numerical analysis, we set cro = 1- The initial guesses used are 
Hc{<j) = ^((T — 0.6) and K{a) = (5(f7 — 0.6) — ^(cr — 0.3). The domain is broken up into bins such that after integration, 
each bin being reassigned the maximum weight. 

We plot the results for Hc{(t) and K{(t) with binsize 0.0025 in Figure [31 Both solutions obey their respective 
constraints thereby justifying the initial ansatz. Using these numerical results we can also compute the average 
microscopic conductivity of the system near the transition. More precisely. 



< a >= ci + C2e 



1/2 



(25) 



with Ci = 0.988 and C2 = 1.945. Note that the maximum conductivity of the system is | in the fully occupied case. 

Monotone sequences: From Figure [31 the domain of Hc{a) appears to bounded away from zero and from the 
maximum value of 8/3. As a check on our numerices, we determine these bounds with the following argument. 
Expanding Hc{a) as a sum of delta functions. 



Q 

Hc{cr)^^A,S{a~a,) < A, < 1, < a, < -, 



(26) 



and iterating Eq. [20lri times yields the sequence, a„,i,a„_2, an,3, •••, which is arranged from largest to smallest. The 
(n + l)th iteration yields the sequence, a„+i^i, a„+i,2, an+i,3, •■•■ 

To determine the largest value of the a^ after the (n + l)th iteration, assume a„.i > a„_„ such that 



1 + Cln.l + an,u 



= 1 

< 1 



1 



1 

2a„ 1 



1 + 2a„,i 



(27) 



We also have 

3a„ 



1 + 3a„ 
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3a„ 


.1 + 1 
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3a„,i n 


Pnniisi 



(28) 

So the largest value of Ui after the {n + l)th iteration is: j^^^^—- Focusing on the largest values of each iteration, we 
form the sequence Q with 

Q = {o-i,i, 0'2,i, 0,3^1, ... a„_i,a„+i^i, ...} 
= {ai,i, /(ai,i), /(/(ai.i)), -, /^""'^ («i,i), /'"^ (ai,i), •••} 

Since an,i > f{an,i) = o.n+1.1 ~ i+a'^ ' Qn.i, sequence Q is a monotonic decreasing sequence with a lower boundary 
of 0. Therefore, sequence Q converges to a finite limit. A, determined by 

with A — 2/3. Therefore, sequence Q converges to 2/3 with the largest value of a^ = 2/3. 
To obtain a lower bound on a^, we construct another sequence, 

Q' = {al,i,g(al,i),...,g("-i)(<i),ff^"n<i),-}, (29) 

where a'^ ^ is the smallest number after each iteration and g{x) = ji^. Since a^ ^ < di'^'n 1) — 1+2'a'^ ' sequence Q' 

is a monotonically increasing sequence with a boundary of unity so that the sequence approaches a finite limit. A' , 
with 

9 A' 

A'^giA') 



2A' + 1 

such that A' = 1/2. Therefore, sequence Q' converges to 1/2, i.e. the smallest value of a^ is 1/2. 

Consequently, the a^s are confined between 1/2 and 2/3 for Hc{<j) as demonstrated in the numerical analysis. One 
can also extend this analysis to Kogut's fc = 3 site percolation analysis to demonstrate that the upper and lower 
bounds on a^ are 1/2 and 1/3 respectively. (In the site formulation of the problem, occupied sites are surrounded 
by half-bonds with some conductivity and unoccupied sites by half-bonds of zero conductivity.) In Figure |4] we plot 
Hc{(j) and K{a) for the k = 3,z — A site percolation problem analyzed by Kogut to demonstrate the bounds. We 
note that the data suggests other gaps in the domains of these functions for the both the bond and site problem. 
Such gaps could indicate a fractal structure. 



III. QUANTUM fc-CORE CONDUCTION 

A. Quantum conduction formulae 

To examine quantum fc-core conduction on a dilute Bethe lattice, we require the conduction formulae for adding 
quantum resistors in series and in parallel. Anderson and collaborators t3Q| have derived the quantum equivalent of 
Ohm's law for two quantum wires in series. Their starting point is the Landauer approach to conductance, g, i.e. a 
scattering matrix approach [311] . Consider two scatterers in series. See Figure [5] The logarithm of the transmission 
probability of the two scatterer system is 



InT = InTi + lnT2 - ln{l + R1R2 - ^/R^cos{0)), (30) 
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FIG. 4: Left: Plot of Hc{cr) with binsize 0.0025 for k = 3,z = 4 classical site percolation. Right: Plot of K{a) for the same 
conditions. 



FIG. 5: Two quantum scatterers in series where fi and ri denote the transmission amplitude and reflection amplitude, respec- 
tively, for the first scatterer, and t2 and r2 the transmission and reflection amplitudes for the second. 



where Ti = |iip with ti representing the transmission amplitude for the first scatterer, etc., and 6 is the phase 
difference between the two scatterers. Assuming 9 is randomly distributed, averaging over the phase difference and 
applying Landauer's formula, we arrive at the conductance formula (in dimensionless quantum units) for 2 scatterers 
in series, 



(1 + -) = (! + -)(! + -), 
9 51 52 



where gi and 52 are typical conductances. In general, for n scatterers in series: 

n -. 



5 



9^ 



(31) 



(32) 



As for the formula for adding quantum resistors in parallel, Arnovas and collaborators [32l ] have demonstrated, via 
a redefinition of a transfer matrix to quantify vertical propagation (as opposed to horizontal propagation in the series 
case), that the transmission probability is interchanged with the reflection probability. Therefore, to 



(1 + 5) = 11(1 + 50 



(33) 



for n quantum resistors in parallel, again, assuming the phase randomizes between scatterers. 

In applying these formulas to the Bethe lattice, we assume one scatterer per occupied bond. Using Eqns. [32] and 
1231 we define the quantum analogs of Sn and T„, or S*^ and T^ as 



n^llil+g,)-!, 



and 



<_>„ 



goTn 



(34) 



(35) 



Tn+ go + l 

where 50 denotes the typical conductance of an individual scatterer. Note that while there is disorder in terms of the 
dilution p, there is also randomness in the individual conductances. The latter disorder is of the Anderson type, while 
the former is of the quantum bond percolation type. 



B. Bounds on pq for fc = 2, ^ = 3 

We begin by assuming that the onset of quantum conduction is driven by the geometric percolation transition. In 
other words, Pq = p^ and 

P''{g)=R5{g) + {l-R)m(g) (36) 

with 

m{g) = 2pRJdg,m{g,)5{g~Sl)+p{l - R) j dg,dg^m{g^)H'i{g^)5{g ^ Si). (37) 

If an occupied bond is connected to infinity geometrically via occupied bonds, it is also connected to infinity quantum 
mechanically via occupied bonds. Assuming the quantum conduction transition is continuous just as the geometric 
transition, we propose 

ff^Cg) = ^H'^{^)- (38) 

Inserting this ansatz into Eq. |37]leads to the e independent terms: 'H'^{s) and Jdxi//^(a;i) exp(— s(a;icro)/(l +o'o))- 
Since these expressions cannot be equated, the initial assumption of the quantum conduction transition being driven 
by the geometry is incorrect, as expected. Loosely speaking, quantum interference prevents extended states in narrow 
channels. In other words, there are not enough occupied bonds participating in the spanning k — 2 cluster at the 
transition to warrant a quantum transition. Certainly, the obvious lower bound on pq is Pc- To compute a better 
lower bound on Pq, one needs to calculate the quantum mechanical version of R with the two types of disorder. 

We now analyze the fully occupied case to determine ii Pq < 1, i.e. if there exists a quantum conduction transition. 
For the fully occupied Bethe lattice that there is a critical value of go below which there is no quantum conduction 
even for p ~ I. This critical value goc is given by 

goJ(l + 5|)^~l] 
l+5oc + [(l + 5|)2-l]- 

For 5| << 1 near the transition (assuming it is continuous), goc = 1/2. For general z, goc — l/(-2 — 2). This result 
agrees with Shapiro [33| . 

We perturb about p = 1. For z — 3 and p — 1, g^c = 1, we, therefore, choose g^ = 2 as an example and eventually 
invoke the expansion parameter c= \ — p « 1. The conductivity g^ of a perfect branch is given by: 

! + - = (! + J)(l + 77— T^)' (40) 

9b 2 [l+gbY-l 

yielding t/f, = 1. Following Stinchcombe [23, we denote g'"^ and g)^'' as the conductances of the branching network and 
of any one of the z branches incident on the origin where one bond has been removed from the nth shell, respectively. 
Note that removing one bond from the nth shell does not initiate the removal of other occupied bonds for fc = 2, z = 3. 
For n>2 

l+g(") = (l + g,)2(l+5(")) (41) 

and 

^ + 4) ^ (1 + ^)(1 + 7 ,, ^ („-i), )• (42) 

Qb 2 (i+<?b)(i+<?r 0-1 



SI = , '':''\ir'.i. .. ■ (39) 



These two equations result in 



1 



3H-3 2 ' g("-i)-f (^^^ 

Starting with an initial value of g^^' — 3, this sequence converges to a finite value of approximately 7. The number of 
ways of removing a bond from the nth shell In > 1) is z{z — 1)"~^. Therefore, the average conductivity with a small 
concentration c of absent conductors is given by 

g{c) = g{c = 0)[1 - c^ z(z - l)«-i(^i^_l_^)]. (44) 

The sum in the above equation diverges as n — >■ cx). However, one can always choose a small enough number c such 
that c^^-^2;(z — ^)'^~^ { ^ ~JZo) — ) remains small. Consequently, there is quantum conduction just below p = 1, 
provided go is large enough. Therefore, pc < Pq < 1 for go > 5oc- 



10 



k = 3,z — 4: case 



We, again, assume that the onset of quantum conduction is driven by the geometric transition. This assumption is 
more plausible for this particular value of k since the onset of the infinite cluster is discontinuous such that one may 
expect a quantum transition. More specifically, we assume (1) Pc = Pq and (2) 



PHg) = RS{g) + il-R)m{g), 



(45) 



where H'^(g) is dictated by the quantum version of Eq. [181 Assuming H'^{g) has the same scaling form as in the 
classical k = 3, z = 4 case, we find Xq — 1/2. We can also construct the lower bound and upper bounds for the domain 
of H^{g) using the same arguments in Sec. II and find an upper bound of (v5 — l)/2 and zero as the lower bound. 
In Figure [6l using the same algorithm as in the classical case to solve the nonlinear integral equations, we plot H^{g) 
and K'^{g). We find that the average microscopic conductance is 



< g >= di + ^26^/^ 



(46) 



with di = 1.379 and ^2 = 4.959. For the Bethe lattice, the microscopic quantum conductivity, < a'^ >, is proportional 
to microscopic quantum conductance and so < cr' > has the same scaling with e. 

The significance of this result is two-fold; (1) we have found a quantum conduction transition with pc — Pq and (2) 
we have found a discontinuous onset of quantum conduction in an Anderson model with two types of disorder, one 
correlated — fc-core correlated — and one not. Therefore, we have discovered a new universality class for an Anderson- 
type transition with the classical random first-order phase transition in the geometry driving the quantum transition. 
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FIG. 6: Left: Plot of H^ig) with binsize 0.005 for k = 3,z = 4 quantum bond percolation. Right: Plot of K''{g) for the same 
conditions. 



IV. DISCUSSION 



In the world of metal-insulator transitions, there are typically two effects to consider: that of electron-electron 
interactions and that of disorder. Typically, disorder-driven (Anderson) transitions [34| . in the absence of interactions, 
are continuous, while interaction-driven transitions in pure systems, such as the Mott-Hubbard transition [35|, are 
discontinuous. Our finding blurs this conventional wisdom in that we have discovered a discontinuous onset in quantum 
conduction as a function of the fc-core correlated disorder in the absence of electron-electron interactions. Presumably, 
there are other geometrical correlations to be constructed and studied — ones that will affect the usually continuous 
nature of the Anderson transition. 

Can this fc-core disorder be realized in an actual experiment? One motivation for fc-core (bootstrap) percolation 
is to ca ptu re some aspect of the principle of local mechanical stability in a static, amorphous packing of jammed 
spheres [22l | . Perhaps a quantum analog of this can be realized in low-temperature packings of metallic nanoparticles? 
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An experiment has already been conducted with a collection of silver quantum dots sitting atop of a Langmuir 
monolayer at room temperature [36|. As the interparticle spacing decreases by compressing the floating particles 
together, the electronic transport goes from hopping to tunnelling to ordinary metallic transport. The authors 
claimed that disorder in the particle size and in the charging energy probably does not drive the transition and, 
instead, argue for a possible first-order Mott transition at room temperature. However, in light of the analysis of the 
onset of classical conduction for A; = 3, we argue for a possible classical correlated percolation transition in conduction. 
Finally, our system contains two types of disorder-the fc-core dilution disorder and the disorder in the individual 
conductances. It would be interesting to retain only the fc-core dilution disorder to realize a fc-core version of quantum 



bond percolation. Perhaps then we would also find Pc — Pq'? We note that Luck and Ashavi [37| have analytically 
investigated quantum conduction on the fully occupied Bethe lattice in the absence of site disorder. They found a 
band structure in the fully occupied case that shrinks to zero as z is increased. They did not analyze the dilute case, 
however. Harris has analyzed quantum uncorrelated bond percolation as the localization transition is approached 
from below by solving the Schrodinger equation ioi E = eigenstates on finite clusters [13, [l3- Since no finite 
clusters exist for fc > 3 on the Bethe lattice, a nontrivial extension of this analysis is required. Such a task, however, 
should be pursued to discover other models with correlated disorder exhibiting novel metal-insulator transitions. 

JMS would like to acknowledge funding from NSF-DMR-0645373 and the hospitality of the Aspen Center of 
Physics where some of this work was completed. The authors would like to thank Ron Maimon for a critical reading 
of an earlier version of the manuscript. 
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